Three-body bound states in atomic mixtures with resonant p-wave interaction 
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We employ the Born-Oppenheimer approximation to find the effective potential in a three-body 
system consisting of a light particle and two heavy ones when the heavy-light short-range interaction 
potential has a resonance corresponding to a non-zero orbital angular momentum. In the case of an 
exact resonance in the p-wave scattering amplitude, the effective potential is attractive and long- 
range, namely it decreases as the third power of the inter-atomic distance. Moreover, we show that 
the range and power of the potential, as well as the number of bound states are determined by the 
mass ratio of the particles and the parameters of the heavy-light short-range potential. 
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Introduction. One of the most intriguingphenomenon of 
few-body physics is the Efimov effect [l|, which mani- 
fests itself in an infinite number of weakly bound three- 
body states if at least two of the three two-body subsys- 
tems exhibit a single weakly s-wave bound state or res- 
onance. The underlying effective potential is attractive 
and decreases as the second power of the inter-atomic dis- 
tance In this Letter, we consider a three-body system 
consisting of a light particle and two heavy ones, when 
the heavy-light short-range interaction potential has a 
weakly bound or quasi-bound, i.e., resonant state with 
a non-zero orbital angular momentum. We show that 
in the case of the exact p-wave resonance, the effective 
potential is also attractive and of long-range, but now 
decreases as the third power of the inter-atomic distance. 

The Efimov effect occurs in systems where the two- 
body scattering length ao is large compared to the char- 
acteristic radius R of the two-body interaction and dis- 
plays an universal behavior, that is the details of the 
underlying short-range physics become irrelevant. In the 
resonant limit, i.e., \oq\ — ¥ oo, the energies E n of the 
three-body states form a geometric sequence, with the 
common ratio determined by the exponent exp(27r/so), 
that is \n{E n /E n+ i) = (2tt/s ) with n = 1,2,3,... The 
parameter so depends on the masses of the particles and 
the number of participating resonant two-body interac- 
tions (two or three). 

Examples of systems with a large scattering length are 
halo nuclei @, 0] and the helium trimer 0] ■ In both cases 
ao is exceedingly large, but not tunable. However, in or- 
der to observe the Efimov spectrum, it is crucial to be 
able to tune ao- In the domain of ultracold atomic gases 
this task is achieved by Feshbach resonances [f| and dif- 
ferent features of the three-body recombination process 
as well as the scattering of the atom off the shallow dimer, 
which are associated with the Efimov effect, have been 
measured Q in this way. An additional prerequisite to 
detect many Efimov states is the use of an atomic mixture 
with heavy atoms of mass M and light ones of mass 
to, since in this case the ratio of two nearest bound-state 



energies is \E n+ i/E n \ ~ 1 for M/m — > oo [§-flfj||. 

In our Letter we consider such an atomic mixture. 
However, in contrast to the standard Efimov scenario we 
focus on an exact p-wave resonance in the heavy-light 
short-range potential and determine the effective inter- 
action potential between the two heavy atoms. More- 
over, we demonstrate that the spectrum of bound states 
is solely determined by the mass ratio of the heavy and 
light particles, and the width of the p-wave state. 
Born-Oppenheimer approach. Our three-body system 
consists of a light particle which interacts with two heavy 
particles and therefore can be easily analyzed within the 
Born-Oppenheimer approximation For this rea- 

son, the Schrodinger equation for the full wave function 
<i>(r, R) = ^>(r; R)x(R) separates into two equations and 
the one for the light particle reads 
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with r± = r ± |R. Here [i = 2mM/(2M + to) « to and 
R denote the reduced mass and the separation between 
the two heavy particles, respectively. For the sake of 
simplicity we assume the potential U to be spherically 
symmetric, i.e., U(r) = U(r), and to have the finite range 
R Q , i.e., U(r > R ) = 0. 
The bound-state energy 



V(R) 
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(2) 



of the light particle serves as an interaction potential for 
the relative motion of the heavy particles given by 
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V(R)] x(R) = 0, 



(3) 



with E being the total three-body energy. 
Interaction potential from scattering approach. Next we 
determine V by a self-consistent scattering of the light 
particle off the two heavy ones. For this purpose, we cast 
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Eq. ([T]) into the integral equation 12 1 



coefficients 
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dr' [U{t'_) + U(t' + )] *(r')' 
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Since the total heavy-light potential i/(r_) + U(r+) is 
nonzero only inside two spheres of radius Rq centered at 
r = ±gR, we represent Eq. @ as the superposition 

*(r) = # (_) (r) + * (+) (r) (5) 
of the two waves 

*<±>(r)= J dr'a^Hr')^^- (6) 
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with 
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into the modified spherical Bessel functions Ti(z) = 
V^7(2l)/ i+1/2 (z) and /Q(z) = y/^J{2z) K l+1/2 {z) (H, 
which is valid for r > r', transforms Eq. ([6j) into 

00 

* (±) W-E E C^M^WnrJ. (9) 

1=0 \mi\<l 

Here Yi mi (n r ) = Y; mi (0 r , <p r ) are the spherical harmonics 
with n r = r/r = (6 r , (p r ). 

We regard the coefficients determined by the in- 
tegral in Eq. ([5]) as independent variables and apply scat- 
tering theory to obtain from Eq. ([9]) explicit equations 
for coupled by the 5-matrix elements of the poten- 
tial U . For this purpose we consider a vicinity of the first 
potential well, that is r = — ^R + x with |x| « Rq, where 
the total solution 



*(-f +X ) = E E *W"I*I)W»*) (10) 
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given by Eq. ([5]) can be expanded into the spherical 
harmonics. Here the radial wave function 

R lmi (K\x\) = C«/C / ( K |x|) + 

OO 

^-1)^1*1) E4™ I} 4i (11) 

i'=0 

is determined by the sum of the two contributions result- 
ing from \&( ± )(— iR + x) defined by Eq. ©, and the 



L=0 
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(12) 



originate from the re-expansion 14] of /Cf(re|x — 
R|)iimj (dx-r) into Yz mi (n x ) with the Clebsch-Gordan 
coefficients cI^lo- 

In order to derive an equation for we cast the 

radial wave Ri mi given by Eq. (1111) into the superposition 

iWK|x|) = Ol(K)/l, (1) (iK|x|) +b l {K)h<i 2) (iK\x\) (13) 

of outgoing and incoming radial waves h\ and with 
amplitudes 

,w^t^4e^( k < (14) 



and 



w = ^E<™'^)4i- 



I 00 



(15) 



l'=0 



The spherical Bessel functions of the third kind and 
m are determined (l3j in terms of /C; and I; as K,i(z) — 
-(iri l /2)h i l 1 \iz) and T,(z) = + hf ] (iz)]/(2i l ). 

Since the amplitudes a; and of the outgoing and 
incoming waves are coupled (ill UBI by the ^-matrix el- 
ements 5/ of the scattering potential U, that is 



ai(n) = Si(iK)bi(n), 



(16) 



we arrive at 



eff + - 1] E ^M)^™, = °- (1 7 ) 



v=o 



Similarly we obtain from the second potential well, cen- 
tered at r = iR, the relation 



ct\ + m*) - 1] £(-i) w '^ m,) M)c, 



(+) 

l' mi 



l'=0 



(18) 

Equations (|17[) and (fT8|) constitute a system of linear 
algebraic equations for determining via Eq. ^ the 
waves i if^ ± \ Its solution is nonzero only if the corre- 
sponding determinant vanishes which provides us with a 
transcendental equation for re — re(i?), and thus for the 
interaction potential V defined by Eq. ([2]). The coeffi- 
cients of these equations are determined by the S'-matrix 
elements of the interaction potential U between the heavy 
and the light atoms. 
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Zero-range limit. In order to test our method, we first 
consider a zero-range potential, for which only s-wave 



scattering occurs and the S-matrix elements read [151 , [I 



2 k 



1/dQ 
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In this case, the system Eqs. (|17l) and (JT5J) reduces to two 



(±) 



algebraic equations for C 00 
only if 



and has non-trivial solutions 



[5o(i/c)-l] %$(kR) = ±1, 



(20) 



with QI$(kR) = [l/{2nR)]e- KR defined by Eq. (TT2 
This condition translates into equation 



1 



i - &op 

for £ = kR with the parameters 

i?0 , R 

a = — and p = — , 
ao Ho 



(21) 



(22) 



and coincides with the equation for the bound-state en- 
ergy obtained in Refs. 0, EH for the case of the 
zero-range potential. 

In the case of a s-wave resonance, that is ao = 0, Eq. 
(|2"Tj) has a solution £ = £* w 0.57 only for the plus sign 
on the right-hand side, which translates into the familiar 
Efimov potential 



_ N+(i?)] 2 
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decaying with the second power of R. 
P-wave resonance. Next we focus on the low-energy 
limit, that is on \E\ <C h 2 /(pR 2 ), or kRq -c 1. In this 
case the S'-matrix elements in each partial wave with the 



orbital angular momentum I can be presented [18 
the form of the effective-range expansion 



It 



5j(»k)-1 = (-1)' 



2k 21+1 



l/ai - {n/2)K 2 + (-l) l + 1 K 21 
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(24) 



The resonance regime in the Z-th partial wave is reached 
when the absolute value of the effective scattering length 



a{\ 3> Rq 1+1 . The effective range n is positive for I > 



and linked [HHH 21 1 to the width of the resonance. 

We now consider the case of a resonant p-wave, that 
is a partial wave with 1 = 1, and substitute the S'-matrix 
elements given by Eq. (JM|) into Eqs. <jTTJ) and (ITU . 
Since So(«k) — 1 ~ Si(in) — 1 ~ kRq and Sj>i(zk) — 1 ~ 
(kRq) 21+1 , we neglect the small terms with I = 2,3, .. and 
arrive at two separate systems of equations with respect 
to the projection m; of the angular momentum. 

For nil — ±1) Eqs. (fTTf and (fT5|) simplify to two equa- 
tions for and have a non-trivial solution only if 



[Si(m)-1] 2l^ 1] (Ki?) =±1 



(25) 



A v ± (p;a,) 




v + (p;a, =0) 
v = (p;a, >0) 
i' + (p;a, <0) 
v""( P ) 



Figure 1: (Color online) Effective potentials V± (p; ai) = 
— [h 2 / (2pRo)]v±(p; ot\) between two heavy particles induced 
by p-wave resonant scattering of a light particle on both sides 
of the corresponding Feshbach resonance. Here v± = (£±/p) 2 
are determined by the solutions £±(p;ai) of Eq. (|26|l for 
\ol\ | = 10~ 3 , /3 = 0.5, and a± = 0; ei = — (a\ //?) is the bound- 
state energy of the light atom in the potential U. The dashed 
line describes the familiar Efimov potential v^ ' = — (£*/p) 2 
defined by Eq. |g3j). 



MD and 2lii 1} (Ki?) 



Since Si(ift) is given by Eq. 
-|[(1 + KR)/{nRf]e- KR , Eq. 
distances i? > 2Rq reads 



5} for £ = and 



+ 



- - ie 

with the dimensionless parameters 



±1 



Oil 



3ai 



and /3 



6 



(26) 



(27) 



In the resonant case, that is o>i = 0, Eq. (|26|) has a 
solution only for the plus sign on the right-hand side. In 
the limit of < £ < 1, we find £ = (Pp)~^ for p » 1, 
which translates into the potential 



y(i,±i) 



[7i«4 1; (.R;0)] 
2^ 



6 



2fi ni? 3 ' 



(28) 



which is independent of Rq. 

The potentials v£ 1,±1} (ii; e*i) = -[/i7(2pJ? 2 )](£ ± /p) 2 
determined by the solutions £± (p; ai ) of Eq. (|26|) are 
presented in Fig. [T] for |ai| = 10 -3 and j3 — 0.5, and 



ot\ =0. The form of V. 



(i,±i) 



is determined by the sign 



of oi, that is by the sign of a\, Eq. (f2"Tl) . Indeed, for 
ct\ > 0, i.e., in the case of the weakly-bound p-wave 
state in U, Vr' ±:l ' as well as V_' ±:l approach for large 
distances, i? > i? |ai| -1 / 3 , the bound state energy e\ = 
— (ai//3) of the light particle. For short distances, i? < 

i? |ai|- 1 / 3 , V { + ,±V) {R;±\ax\) approach v£ ,:kl) (R; a x = 
0) =V^ ±1 \R). 
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In the case of a p-wave resonance, the matrix element 
Si corresponding to the resonant channel is of the same 
order as So for the non- resonant channel llj. Therefore, 
for mi — we have to take into account in Eqs. (fTT|) and 
(fT5|) both the s- and p- waves, which gives rise to a system 
of four algebraic equations for Cq*q and , leading us 
to the relation 



1 



(s -i)(Si-i) 



T (50-1)2$ 



(0) 



{Si IWn 



(29) 



According to Eqs. (T2J|, dHJ) and J27|, we obtain a o ° } : 
21$ = -a^[(l + K R)/( K R) 2 }e- KR and = f^i? 2 - 



2KR + 2)/{nR) 3 ]e- KR , and Eq. 
the form 

p-2? 



for £ = k^R takes 



(e- ao p)(^ 2 -«iP 3 -ie 3 
(e 2 + 2e + 2) i 



/3,o£ 2 - «ip 3 



if 3 

3 ^ 



£ - "op 



(30) 



In the resonant case, ai = 0, we have ao ~ 1 and Eq. 
(|30| has a solution only for the plus sign on the right- 
hand side. In the limit of < £ -C 1, we find £ = 
[2 /(Bp)]? for p > 1, giving rise to the potential 



(1>0) ^ [^f(i?;0)] 2 

2/z 

The potentials v£' 0) (i?; ai) = 



fi riR 3 



(31) 



-[ft 2 /(2/^)](£ ± /p) 2 
determined by the solutions £±(p;ai) of Eq. ([50)1 at 
ao ~ 1 and /3 ~ 1 are similar to V±' of Fig. [I] with 
two qualitative differences: (i) the ranges of v!£' ±1] and 
Vj. 1,0 " 1 are different and equal to Ri = i?o|cn| -1 ^ 3 and 
i?2 = i?o(0.5 |«i I)" 1 / 3 , respectively, and (ii) in the case 
of exact resonance, V^ 1,0 ) given by Eq. (f5T|) has the same 
asymptotic behavior as y' 1 ' ±1 - ) defined by Eq. ((28)) with 
twice the amplitude. 

Spectrum of induced 1 / i? 3 -potential. Finally we focus on 
the dynamics of the two heavy particles dictated by the 
Schrodinger equation (|3]) with the potential V given by 
Eqs. (|28p and (l3~Tj) and induced by the p-wave resonance 
in the light-heavy interaction. We emphasize that V is 
only meaningful for R ^> i?o , since for R ~ Rq it is likely 
to be determined by the inter-atomic forces and obviously 
cannot be treated by the effective-range expansion. 

The energies E n of the bound states with zero angular 
orbital momentum follow from the familiar WKB quan- 
tization rule 11 1 



n — uq 



^ J y/M[E n - V(R)] dR, 



(32) 



giving rise [22, 23 1 to the spectrum 

ft 2 



MRl V 9 



no — n 



(33) 



for the weakly bound states. Here we have introduced 
the characteristic range 



j^ = (2-kd — 



(34) 



of the effective potential V in the resonant case, that is 
ai = 0, and [no] — n = 1,2,.., where the integer part 
[no], determined by the phase of the wave function at 
the short distances R ~ Ro, plays a role of a three-body 
parameter and g — r(|)/[y / 7rr(|)]. 

Since V defined by Eqs. (|28l) and (l3"Tj) has a tail falling 
off faster than — 1/i? 2 , it supports [ll| only a finite num- 
ber iVo of bound states with zero angular orbital moment. 
Indeed, Nq can be estimated [ll[ by the WKB method 
and yields 



N 




(2-|m/|)3M 



nRo 



m 



(35) 

that is No is determined by the square root of the ratio of 
the mass-ratio M/m to the dimensionless effective range 
riRo of the p-wave resonance. 

Scattering off ' 1 / R 3 -potential and summary. The appear- 
ance of V given by Eqs. (|2"8]l and ([3"Tjl can be verified ex- 
perimentally by scattering a heavy atom off the diatomic 
molecule consisting of a heavy and a light atom. The 
predicted three-body bound states manifest themselves 
as resonances in the cross-section of the atom-molecule 
scattering when we tune the magnetic field close to the p- 
wave Feshbach resonance. Moreover, due to the inverse- 
cube tail the cross-section Oh of the L-th partial wave has 
the unique behavior a (E) = irRl\n(MRlE/h 2 ) 
and o- L >o(E) = 7ri? 2 [(2L + 1)/(L 2 + L) 2 ] at the low inci- 
dent energy E < h 2 /{MR 2 ). 

In summary, we have found a novel series of bound 
states in the three-body system consisting of a light parti- 
cle and two heavy ones when the heavy-light short-range 
interaction potential has the p-wave resonance. In the 
case of an exact resonance, the effective potential is at- 
tractive and of long-range. Moreover, the spectrum of 
bound states is determined by the mass ratio of the heavy 
and light particles as well as the parameters of the heavy- 
light short-range potential. 
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